Screening Behavior and Scaling Exponents from Quantum Hall Wavefunctions 
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We provide a robust and generic method to assess the screening properties and extract the scaling 
exponents of quasiparticle edge excitations of quantum Hall states from model wavefunctions. We 
numerically implement this method for the fundamental quasihole and hole excitations of several 
model states. For the Laughlin, Moore-Read, and Z3-Read-Rezayi states, we find agreement with 
the predicted edge theory, verifying the bulk-edge correspondence. We also use this to obtain the 
first clear microscopic demonstration of the pathologies of the Gaffnian wavefunctions. 
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Since Laughlin's explanation [T] of the fractional quan- 
tum Hall (FQH) effect ;2 in terms of incompressible 
quantum liquids, model (trial) wavefunctions have played 
a significant role in understanding FQH states. This 
understanding has benefited greatly from the realiza- 
tion J3] 2] that conformal field theories (CFTs) [3] can 
generate model FQH wavefunctions (including Laugh- 
lin's) using conformal blocks of appropriately chosen 
CFT operators. The analytic properties of such wave- 
functions are manifested by the CFTs from which they 
are produced, allowing straightforward extraction of im- 
portant experimentally measurable quantities character- 
istic of the universality classes which these wavefunctions 
(are intended to) represent. Such quantities include the 
charges and braiding statistics of quasiparticles, as well 
as the edge modes' scaling exponents and chiral cen- 
tral charge, which are measurable, for example, through 
tunneling, interferometry, and thermal transport experi- 
ments [6HS]. 

However, the validity of a CFT based FQH state 
and its experimental predictions depends upon the as- 
sumption that it produces legitimate FQH wavefunc- 
tions whose inner products (with respect to particle co- 
ordinates) match those of the CFT. This assumption 
has been verified for the Laughlin pQ and Moore-Read 
(MR) [4] states by mapping their respective wavefunc- 
tions' inner products to the partition functions of certain 
classical two-dimensional plasmas [TJ [§] that are in their 
metallic phases [TO] [IT] . More generally, plasma map- 
pings can be constructed for proposed FQH states based 
on all the minimal model CFTs [9] and their general- 
izations [12] , including the Read-Rezayi (RR) states [13] , 
but the corresponding plasmas are more complicated 
and their screening properties not yet well-established. 
Hence, the question of which CFTs produce legitimate 
FQH wavefunctions remains largely unsettled. 

One may also consider states generated using non- 
unitary CFTs, such as the "Gaffnian" (Gf ) state [H IT5] 
generated from the M(5,3) non-unitary minimal model 
CFT [together with a U(l) Laughlin-like charge sector]. 



It is difficult to envision how FQH wavefunctions could 
emerge from non-unitary CFTs free of pathologies, such 
as non-unitary edge modes and braiding statistics or gap- 
lessness in the bulk, and arguments have been made |15l - 
[T7] that they cannot describe gapped phases. However, 
the failures of such states can be subtle and a com- 
plete microscopic understanding of their deficiencies is 
still missing. Numerical studies of Gf wavefunctions for 
small system sizes produce some indications that the Gf 
could be an acceptable FQH state [TS1 [TM2T)] , while the 
first signs of problems appear in numerical studies of the 
entanglement spectra |21j . which suggest non-unitarity 
in the edge modes, and of states with quasiholes |18) . 
which exhibit indications of level crossings when deform- 
ing between the Coulomb and Gf Hamiltonians. How- 
ever, these results give no clear indications of the behav- 
ior in the thermodynamic limit and clear demonstrations 
of the anticipated failures of the Gf wavefunctions have 
not been identified. 

In this paper, we provide methods of analyzing model 
wavefunctions of CFT generated states that can help 
determine whether they exhibit the required properties 
of FQH states and verify the bulk-edge correspondence, 
which conjectures that the (1 + 1)D edge theory is de- 
scribed by the same CFT used to generate the 2D bulk 
wavefunctions, from microscopic computations. We de- 
vise a numerical method of extracting the scaling expo- 
nents of edge excitations from microscopic wavefunctions, 
using Jack polynomial expansions to produce sequences 
that must converge in the thermodynamic limit to the 
exponents for a properly screening state. Applying these 
methods, we find behavior consistent with the respective 
CFT descriptions for the Laughlin, MR, and RR states 
(though system size limitations prevent a clean extrac- 
tion of the RR state's quasihole exponent). We analyze 
only bosonic FQH states here, but the generalization to 
fermionic states is straightforward and yields similar re- 
sults [22] . We find the behavior for the Gf wavefunctions 
to be inconsistent with proper screening and the CFT 
description of a FQH state, providing clear indications of 
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their failure. 

The basic configuration that we use to study candidate 
states is the planar disk geometry with some number of 
quasiparticles, one of which is located at (complex) co- 
ordinate i], while all others are located at the origin of 
the disk. The corresponding (unnormalized) wavefunc- 
tion for N bosonic particles with coordinates z\ , . . . , zn 
can generally be written in the form 
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where P a (zi, . . . , zm) are symmetric polynomials, is 
the highest power of rj, and I — \JhcjeB is the magnetic 
length. The inner product of the model wavefunction 
with its quasiparticle at possibly different positions is 



/N nj, 
.7=1 a=0 



J = l 

where we have defined 
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to be the norm-squared of P a (zi, . . . , Zat), and used the 
fact that these polynomials with different a are orthogo- 
nal, since their total orders are unequal. 

While a quasiparticle cannot physically exist outside of 
the quantum Hall droplet, one can formally write model 
wavefunctions with the quasiparticle coordinates located 
anywhere. The configuration with a quasiparticle out- 
side the Hall droplet can be used to assess the candidate 
state's screening properties, extract the quasiparticle's 
scaling exponent, and determine whether the bulk-edge 
correspondence holds, using arguments similar to those 
applied to Laughlin states by Wen [5] . If the state screens 
properly, the norm of the wavefunction with the quasi- 
particle outside the disk should take the form 



\\*f = mv)~C\r)\ 2n * 
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in the thermodynamic limit (N —> oo) for \r)\ — R £, 
where C is an overall constant, R is the radius of the 
disk of Hall fluid, and the exponent g is a constant that 
depends on the state and quasiparticle. This form can 
be justified for certain model states (including the ones 
studied here) by mapping the norm to the free energy of 
a classical 2D plasma [U [9l [22] . In its screening phase, 
the plasma behave like a metal, so the free energy of the 
system is the Coulombic energy of a metallic disk with 
a point charge outside it. This can be computed by the 
method of images, which requires an image charge at 



R? jr\ that gives rise to the important f 1 — iff 7 ) term 
in Eq. Q. Analytic continuation of Eq. Q gives 



T(fj,r]') ~ CWf'E 
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where the approximation holds in the regime where 
\v\ i \v'\ ^> R, and the thermodynamic limit gives an infi- 
nite sum that converges for |r/| , \rf\ > R. 

If the candidate state is properly screening and gives 
rise to a well-defined edge theory described by a CFT, 
then the inner product of states with a quasihole on the 
boundary of the Hall droplet should match the (equal 
time) quasihole correlator in the edge CFT. Specifically, 
for -q = Re %e and rf = Re 10 in the 2D planar wavefunc- 
and w' = eTi +iB ' in the (1 + 1)D 



tion, and w 

edge theory, one should have (up to overall phases) 



T(fj, rf) - ($ f K) $ H) = (to' - to) 
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where h is the conformal weight of the quasihole oper- 
ator $ in the edge CFT. Thus, if these properties hold 
for the candidate state, then g in Eq. ^ should equal 
the scaling exponent g = 2h of the quasiparticle, in the 
thermodynamic limit. 

By matching powers of f)r\' in Eqs. Q and jHj), we find 
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relating the symmetric polynomials' norms to g. Dividing 
by the n = expression (which is just C — Af n<t> ), we 
obtain the sequence of equations 
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defining , an approximation of g coming from the 
(n^ — n)th order term of the polynomial expansion. This 
provides a robust method to numerically compute g and 
provides a consistency check on the properties of the can- 
didate state, since each must independently converge 
to the same value in the thermodynamic limit. Deal- 
ing with finite system sizes, the most relevant expansion 
terms outside the disk are those with n small and the 
accuracy of the approximations will decrease for larger 
n. [The terms in Eq. ([!]) with n large must accurately 
describe the regime where the quasiparticle is inside the 
disk, which exhibits different behavior.] Hence, we focus 
on the n = 1 and 2 expressions 

■A/^-i 
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The model wavefunctions of the states we focus on 
can be simply expressed in terms of Jack polynomi- 
als |23j . which makes them easier to generate numeri- 
cally [21] • The (k,m) series of Jack states at filling fac- 
tor v = k/m includes the Laughlin = (1, 2), MR = (2, 2), 
Zfe-RR = [k, 2), and Gf = (2,3) states. The Jack states 
are comprised of symmetric polynomials that vanish as 
the mth power of the separation of k + 1 particles ap- 
proaching each other. 

We consider a (k, m) Jack state with N particles 
(where N is a multiple of k), a flux quasiparticle 
at the origin, and a fundamental (flux 1/fc) quasiholc at 
position rj. (One can alternatively view the quasiparti- 
cle at the origin as a composite of k — 1 fundamental 
quasiholes.) In this case [19], = N/k and 



P a (zi, . . . , z N ) = (-k) " a J* a (z 1 ,...,z N ) 
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where J" denotes the N particle Jack polynomial with 
Jack parameter a = — (fc+l)/(m— 1) for the (A:,m) series 
and root configuration (in terms of occupation numbers) 



pb a = (l,fc-l,0 



m-2\ a 
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(13) 



e.g. /j = [0,M m ~\fc,0 m -\£,- 
For this configuration, we use R 



FIG. 1: The fundamental quasihole scaling exponent versus 
the number of particles for the model states: (a) v = 1/2 
Laughlin, (b) v = 1 MR, (c) v = 3/2 Z 3 -RR state, and (d) 
v — 2/3 Gf. and g^ values are denoted using squares and 
triangles, respectively. The dashed lines indicate the values 
Sqh = | , | , jq , and predicted from CFT for the respective 
states. The solid lines are quadratic fits in 1/N, which give 

(n) 



, 0'" 1 , fc, m 1 ,k] . the N — > oo extrapolated values, g^ 
= ^J2{v- 1 N + I)l 



for the radius of the Hall droplet, as defined by the 
(neutralizing background) disk of uniform charge den- 
sity p — ve/2Ttt 2 whose total charge is (N + v)e, equal 
to that of all quasiparticles minus all particles. (Note: 
There is some ambiguity in how to define the radius of 
the Hall droplet, as it cannot be a sharply defined quan- 
tity. Alternative choices produce similar results, as their 
differences vanish in the thermodynamic limit.) 

The v = 1/m Laughlin state has the property 
Afn^-i — 2-/W 2 A< r n</) , which can be proven using J^ n _ x — 

Y,%=i z i J Z n<t> -> integration by parts, and Yh=i 557 = 

0. This provides an exact expression = (m + l/N)^ 1 , 
which rapidly converges to the expected value g q h = 1 jm 
for Laughlin quasiholes. 

We numerically compute and gty for the Laugh- 
lin (up to N = 16), MR (up to N = 22), Z 3 -RR (up 
to N — 27), and Gf (up to N — 18) states, as shown 
in Fig. [l] These computations involve squeezed Hilbcrt 
spaces as large as 5.3 x 10 9 . Motivated by the exact result 
for the Laughlin state, we use quadratic fits in 1/N to 
the g^ points with N > 5k in order to extrapolate the 

~(n) 

data to the N — > oo limit. The extrapolated values, , 
are indicated (with fitting errors) in the figure. For the 
Laughlin and MR states, gQ and g^ converge towards 
each other and produce excellent agreement with the pre- 
dicted exponent values (within 2%). As k increases, the 
behavior becomes worse. For the RR state, g$ and g^ 
converge towards each other, but produce gj$ that differ 



from the predicted value by around 20%. We conjecture 
that this discrepancy is due to significant finite-size ef- 
fects arising from the large size of the quasiparticle at the 
origin and that it will vanish as system size increases. In 
contrast to these states, g^ and gfi for the Gf state 
are very far from the predicted value and appear to be 
diverging from it, as well as from each other. 

We note that Hu et al. 25 proposed a different method 
to extract scaling exponents from quasiparticle tunneling 
between edges of an annulus. They found excellent agree- 
ment with the predicted values for the Laughlin-type 
(flux 1) quasiholes of the Laughlin, MR, and RR states. 
However, their extracted exponents for the (non-Abelian) 
fundamental quasiholes of the MR and RR states exhibit 
a (possibly systematic, but so far unexplained) large dis- 
agreement with the predicted values. 

We next consider a (k, m) Jack state with N particles 
(where 7V + 1 is a multiple of k), with no quasiparticles at 
the origin, and a hole (i.e. a flux v quasiparticle with 
the same statistics as the particles) at position r\. (One 
can alternatively view this state as the N + 1 particle 
ground state with one of its particles pinned at posi- 
tion t), i.e. V(t];zi,...,zn) = *o(??, z\, ■ ■ ■ , ^jv)e |r ' |2/4f2 .) 
In this case, n = f {N + 1) - m, R = y/2v~^ (N + 1)£, 
and P a is generally equal to a sum of Jack polynomials, 
rather than a single Jack. 

To obtain the polynomials P„ . , P n . _i, P n , _2 needed 
to compute £t an d g^ f° r l ar S e sizes we require more 
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FIG. 2: The hole (and particle) scaling exponent versus the 
number of (unpinned) particles for the model states: (a) v = 
1/2 Laughlin, (b) v = 1 MR, (c) v = 3/2 Z 3 -RR state, and (d) 

v = 2/3 Gf. an d 5h 2 ' values are denoted using squares and 
triangles, respectively. The dashed lines indicate the values 
gh = 2 predicted from the CFT for the unitary states. (CFT 
predicts <?h = 3 for Gf.) The solid lines are quadratic fits in 
1/N, which give the N — > oo extrapolated values, g^ 1 . 



refined Jack machinery and simply quote the results. In 
all cases, we have 
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For the Laughlin, MR, Z3-RR, and Gf states, the coef- 
ficients in are given by a\ — 1, 1, 1, and —3/2, 
respectively, and a 2 = 14/5, 5/7, 26/81, and 21/16, re- 
spectively. 

Similar to the quasihole case, the v = 1/m Laughlin 
state has the property J\f n ,_i = 2N£ 2 Af n<l> , which pro- 
vides an exact expression — m(l + 1/N) -1 , which 
rapidly converges to the expected value <?h = m for 
Laughlin holes/particles. 

We numerically compute g^' and g^ for the Laugh- 
lin (up to N = 16), MR (up to N = 21), Z 3 -RR (up 
to N — 26), and Gf (up to N = 19) states, as shown 
in Fig. [2j These computations involve squeezed Hilbert 
spaces as large as 1.5 x 10 10 . We use quadratic fits in 1/N 



to the g^ points with N > 5k — 1 in order to extrapolate 
the data to the N —> 00 limit. The extrapolated values, 
ffj; , are indicated in the figure. For the Laughlin, MR, 

and RR states, we see that g^ and g^ converge toward 
each other and toward the value g^ = 2 predicted from 
CFT, producing excellent agreement between g^ and 
the predicted value. For the RR state, these are signifi- 
cantly better (within 6% of the predicted value) than the 
results for quasihole exponents, likely because of the ab- 
sence of (large) quasiparticles at the origin. Surprisingly, 
the Gf state also exhibits nice convergence of g^ and g^ 
toward each other and toward a value g^ « 2.37, which 
differs from the value g^ = 3 predicted from CFT. More- 
over, it violates "spin-statistics," which requires bosons 
to have even integer values of g (and non-unitary spin- 
statistics, which would allow odd integer values of g for 
bosons) . 

We have provided a robust method of testing the 
screening properties of candidate FQH states, extract- 
ing the edge excitations' scaling exponents, and verifying 
the bulk-edge correspondence from microscopic model 
wavefunctions. Applying these methods for fundamen- 
tal quasiholes and holes, we find that the Laughlin and 
MR states behave as they should for a properly screen- 
ing FQH state, with scaling exponents matching those 
predicted for the expected edge CFT. The Z 3 -RR state 
exhibits stronger finite size effects for the quasihole ex- 
ponent, though it appears to have proper screening and 
the hole exponent matches the CFT predicted value. In 
contrast, the quasihole and hole exponent computations 
provide a clear indication that the Gf wavefunctions do 
not exhibit the proper screening necessary for a legiti- 
mate FQH candidate. Despite this, the particle exponent 
computation indicates that there is some sort of screen- 
ing occurring in the Gf state, which may explain why 
the ground state appeared to exhibit qualities similar to 
a FQH state, despite the fact that it, nonetheless, pos- 
sesses pathologies. 
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